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Abstract 

We introduce an extension of the ^-reduced KP hierarchy, which we cah the l- 
Bogoyavlensky hierarchy. Bogoyavlensky's 2 + 1-dimensionaI extension of the KdV 
equation is the lowest equation of the hierarchy in case of £ = 2. We present a 
group-theoretic characterization of this hierarchy on the basis of the 2-toroidal Lie 
algebra st^™. This reproduces essentiahy the same Hirota bilinear equations as those 
recently introduced by Billig and lohara et al. We can further derive these Hirota 
bilinear equation from a Lax formalism of the hierarchy. This Lax formalism also 
enables us to construct a family of special solutions that generalize the breaking 
soliton solutions of Bogoyavlensky. These solutions contain the A-soliton solutions, 
which are usually constructed by use of vertex operators. 

1 Introduction 

In the early 1980s, Date et al. [0 discovered the remarkable fact that the Lie algebra 0[(cx3) 
acts on the solutions of the KP(Kadomtsev-Petviashvili) hierarchy. The link between 
infinite-dimensional Lie algebras and soliton equations has been widely extended. In 
particular, Kac and Wakimoto ||1^ developed a scheme to construct hierarchies from 
vertex operator representations of general afiine Lie algebras. On these backgrounds, 
Billig and lohara et al. derived new hierarchies of Hirota bilinear equations from 



2-toroidal Lie algebra g*""", a central extension of double loop algebra 0[s^^, t^^] of a simple 
Lie algebra q. Applying the construction to the case g*°'' = st"^^ (and their principal vertex 
realization), we have an extension of the ^-reduced KP hierarchy of Hirota bilinear form. 
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For the simphcity of exposition, we here deal with the case g = sl2- The typical Hirota 



bilinear equations of lower degrees in hierarchy of |P1,[T5] are the following 



(D^-4D,A)r-r = 0, (1.1) 
{DyDl + 2DyDt-QD,D,)T-T = Q. (1.2) 

If we count the degrees of Hirota's D-operators Dy, D^, D^, Dt as 0,1,2,3 respectively, 
the Hirota bilinear equations above have degrees 4 and 3 respectively. The first one (|1 . 1|) 

is a famous bilinear form of the KdV equation Aut = 6uUx + Uxxx, where u *== 2(logr)a;a;. 
The second one (|1.2| ) can be written as a non-linear equation known as Bogoyavlensky's 
2 + 1-dimensional equation 

Uz = '^Uxxy + VUx + UUy, (1.3) 

where we set v =^ (logr)^,^^ and used the KdV equation to eliminate the terms including 
derivative dt- Equation ( |1.3| ) has the following equivalent Lax form 

^ = [Pdy + C, P] where P = dl + u,C = vd^ + ^Uy. (1.4) 

The main purpose of this paper, motivated by the fact above, is to establish a Lax 
formalism of the Hirota bilinear equations arising from sll°^ for any i > 2. For this 
purpose, we shall introduce the i-Bogoyavlensky hierarchy. We shall also construct a 
family of special solutions that generalize the breaking solitons studied by Bogoyavlensky. 

We shall develop our discussions as follows. In section H, we discuss the hierarchy on 
a viewpoint of representation theory. First we give definition of the Lie algebra g*™, and 
then explain a general result on their representations(Lemma 0). Using this lemma, we 

construct a representation of sl^ , based on the ^-reduction procedure of 0[(oo) in 0. We 
shall show that every vector r in the SL^^-orbit of the vacuum vector satisfies a hierarchy 
of Hirota bilinear equations, which is essentially the same as those derived in [3,||14||, 



although the generating functions are apparently different. In section |^, we first give a 
review on the Lax formalism of the KP hierarchy. Next we present a heuristic introduction 
to Bogoyavlensky's hierarchy. These arguments lead us to formulate the i- Bogoyavlensky 
hierarchy for each i > 2, where the so-called Wilson-Sato operator plays the significant 
role. In section ^, we give a residue formula for the formal Baker- Akhiezer functions of the 
hierarchy by the calculus of pseudo-differential operators. The formula leads to exactly 
the same system of Hirota bilinear equations derived in section ^ Note that the most part 
of sections can be read independently from the preceding section. Section |^ is devoted 
to investigate the special solutions of the hierarchy, namely the Wronskian solutions and 
the iV-soliton solutions. As an appendix, we give a table of Hirota equations of low degree 
in section ^ 

It should be mentioned that the Bogoyavlensky's equation is a dimensional reduction, 
to 2 + 1 dimensions, of the four dimensional self-dual Yang-Mills equation. Mason and 
Sparling [BH] pointed out that many integrable systems in 2 -f- 1 and lower dimensions 
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(in particular, 1 + 1 dimensional soliton equations such as the KdV equation) can be 
thus derived from the self-dual Yang- Mills equation. In fact, Schiff rediscovered 
Bogoyavlensky's equation, having been unaware of Bogoyavlensky's work. Schiff 's work 



was further extended by Yu et al. pO, [31 



Throughout the paper, the base field is the field C of complex numbers. The symbol Z 
stands for the set of integers. The binomial coefficients (^) is defined as n{n — 1) ■ ■ ■ (n — 
k + l)/k\ for k > Q. For a vector space V, we denote by V[z^\ V[ z^^^w^^] the spaces 

of Laurent polynomials with the coefficients in V 
respectively. For an integer £ > 2, we denote by gl^ (resp. sk) the Lie algebras of the all 
(resp. traceless) i x i matrices. 

Note added: After the first version of this paper had written, we learned that F. 
Calogero had studied equation ( |1.3|) in 1975. We thank Kouichi Toda and his collab- 
orators for the information. 
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2 Deriving hierarchy from representation theory 
2.1 Definitions of the toroidal Lie algebra 

Let be a finite-dimensional simple Lie algebra over C. Let R be the ring of Laurent 
polynomials of two variables C[s^^,i(:^^]. The module of Kahler differentials Qn of R is 
defined with the canonical derivation d : R ^ VLr. As an i?-module, VLr is freely generated 
by ds and dt. Let ~ : VLr — > VL^/dR be the canonical projection. Let /C denote flji/dR. 
Let (•!•) be the normalized Killing form ([1^) on g. We define the Lie algebra structure 

on q'°' = ® © /C by 



[X®f,Y®g] = [X, Y](»fg + {X\Y){df)g, [JC, 0*°^ = 0. (2.1) 



In the paper pO[, Kassel proved that the bracket defines a universal central extension of 
® i? as a Lie algebra over C, for a large class of commutative algebra R. One can find 
in 



24] a simpler proof (a natural extension of Wilson's discussions |^2|) of the Kassel's 
result that works for the cases when the ground field of R has characteristic 0. 

Let V = Rdiogs © Rd\ogt be the Lie algebra of derivations on R. A derivation 5 & V 



can be naturally extended to a derivation on the Lie algebra © i? by 

5 (X © /) = X © 5f. (2.2) 
It is known that a derivation 5 G acting on the Lie algebra © i? has a natural extension 



to 0*°''. The action on the center K, is given explicitly as follows |13|. First we shall define 
the action of T) on VLr by: 

fdiogu{gdlogv) = {fdiogug)dlogv + Su,vgdf for u,v = s,t. (2.3) 

The action preserves the exact forms dR, namely we have fdiogu{dg) = d{fd\ogug) for 
u = s,t, and so induces an action, also denoted by the same notation, on /C = Qji/dR. 

Let V denote the Lie subalgebra Rd\ogt of V. We shall add the derivations V to 0*°'' 
to get the Lie algebra 

0*°- = © V. 

Here the bracket in 0*™ is given by ( ^.11 ). We define the bracket between T) and 0*™ by 
the action (|]2D, and (|]3D, [5, a] = 6{a) {5 et>, ae Q^°'). The bracket in V is defined by 



[fdiogt,gdiogt] = ifidiogtg) - g{diogtf)) diogt - {diogtg)d{diogtf). 

Note that if /, (7 G C[)f:^^] in the above formula, then it is equivalent to the relation of the 
Virasoro algebra. Remark that 0*™ is not the semidirect product of 0*™ with the action 

of:Dby (iD,(p|). 
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We have, for u — s,t, the Lie subalgebras 

with the brackets given by 

[Xu"', yii"] = u^'^'^iX, Y] + mS^+n,o{X\Y) K^, 

which are isomorphic to the affine Lie algebra g with the canonical central element Ku 
dlogu. 

We prepare the generating series of 0*°^ as follows: 

is:^(^) ^Jf ^s»i"^(ilogs- ir^(z) = ^s"t"^cilogi • 

neZ neZ 



where m e Z and A e 0. The relation d{sH'^) — can be neatly expressed by these 
generating series as 



2.2 A Heisenberg subalgebra D of g 



tor 



The Lie algebra 0*°'' has a Heisenberg subalgebra D with basis {s"9iogt, s'^d log dlog s (n e 
Z)}. In fact, if we set 0„ =^ s'^dlogt, 0* =^ s"^ogt, and /i =^ dlogs, then we have 

[0m > 0n] = mSm+nflfl, [<Pm, 0n] = [0m> C] = ^ ] = 0- 

The Heisenberg algebra D is degenerate in the sense that [0q, 0o] = 0. We also note that 
D commute with Qg. 

Next we realize the action of D on the space of polynomials 

= C[ye, y2e, yu, . . . ; yl y^, y^, ...\® C[e^°, e-^°], vac^ = 1 e 
where the action of D is defined by: 



d 



m>0 ^ ^ m > 



^^-^ , ^ <^ ^^-^ ~ , ^ = id. 

-my_rni m<0 (-myl^^ m<0 

For convenience' sake, we adopted rather peculiar notation for the variables yni etc. 
We define the generating series 



and for each m G Z, the following vertex operators 



n>0 




VUz) = exp m 5^ y^,z- e'"^" exp -m -Tpr^"" • (2-4) 



We also introduce the normal product ■.(f){z)Vm{z): '= (j){z)^oVm{z) + Vm{z)(f){z)yo where 

<P{z)<o = En<o and 0(^)>o = E„>o 0^2:"""^ 

We shall say a representation V of the affine Lie algebra g^, is restricted if and only if 
for all u G and X G g, there exists N such that Xs^ ■ v = for Vn > iV. Now we can 
state the following lemma due to lohara et al. and Berman-Billig 



Lemma 1 Let (V, vr) be a restricted representation o/fls such that d\ogs idy. Then 
we can define the representation tt*™ of on V ^ J-'^ such that 

Am{z) ^A^{z)Vm{z), K^{Z) ^ Vm{z) , 

Dl{z) ^c^\z)Vm{z), Kliz) ^:(t>iz)VUz):, 
where A e Q, m G Z and A'"{z) =^ Zlnez 7c{As"')z^"'^^ . 

2.3 Review on a representation of 0[(oo) and its ^-reduction 

In the subsection, we collect basic results on representation of gl(cxD) which is used in the 
sequel. One can find in the original paper M the proofs we omit. We shall mainly follow 



the notation used in book [16]. 



Let the associative C-algebra A be generated by ipi^ifjlii G Z + i) with the relations: 

+ = + = o, ip-^pj + ipjip* = 5^+^- o- 

Consider a left ^-module with a cyclic vector vac^ satisfying 

ijjivac^A = ip*va.c_A = for i > 0. (2.5) 

This ^-module ./lvac_4 is called the fermionic Fock space, which we denote by JF4. 
We define the following infinite formal Laurent series of the variable A 

and the normal ordering of the following quadratic expression to be 

' ' ^ 1 -tp*ilji if i> or j < ' 
Consider the following set of infinite complex matrices 



for \i - i\ > o| . 



It forms a Lie algebra over C with the usual bracket [A, B]q =^ AB — BA, which we shall 
denote by 0[(cx3). We shall define the 2-cocycle u on g[(oo) by : 

u;{Eij,Ei>j>) = 6ji'6if{6{j) - 9{i)) 

where Eij =^ j/g^+i matrix unit and 9 is defined by 9{i) =^ 1 if z > and 

9{i) =^ if 2 < 0. Let us introduce the central extension 0[(oo) =^ sl(c)o) © Cc with the 
central element c and the bracket: 

[A, B] = [A, B]o + u{A, B)c, A,B e ^(cx)). 

Lemma 2 /0/ The map Eij i— >: ip-iip* :,c h-> idjr^ defines the representation tToo of qI{oo) 
on the fermionic Fock space JF4. 

In terms of the generating series defined by 

the lemma has the equivalent expression 

E{X,f^)^:^{X)r{^i): . 
Fix an integer i > 2. Let A be a matrix {aij)ij^^_^_i in 0l(oo). We assume the condition: 

ttij = Qi+ej+e for all i,j eZ + ^. (2.6) 
Then we introduce, for each n G Z, the i x i complex matrix by 

and define an i x i matrix of Laurent polynomial of the variable s 













Conversely, let ^„ ^nS" 


be 


a Laurent polynomial 


in 




-periodic infinite matrix in 


the following block form 










••. ••. ••. • 


. \ 








■•• Ao A ■ 










■•• A_i Ao Ai 




G 0[(oo) 


n 














i^-. ••. ••. ••. • 


■■) 
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Let us denote by t£ the above map, which is clearly injective. Let X,Y G gl^ and consider 
Xs"^,Ys"' {171,11 G Z) as elements of g[(cx3). Then the commutation relation is given by 

= [X,F]s'"+" + m5„+„,otr(XF)c. 

So, in particular, if we extend Li to sle = sle[s^^] © CKg by ie{Ks) =^ c, then it gives the 
embedding li of Lie algebras sli into 0l(oo). 

We shall describe a base of Le^sk) = sli for the later use. 

Definition 1 For each i-th root of unity ( ^ 1, we define G 0t(c>o) (n G Z) by the 
following series 

i?(A,CA) = J]i5;^A— \ 

It is easy to see that satisfies ^-periodic condition (|2.6|). If we shall regard E^ as an 
element of then E^ is traceless. Note that E^ is homogeneous of the principal 

degree n, here by the principal degree of the matrix s^Cij in qI^ we mean j — i + ni & Z. 

Example 1 For £ = 2, ( = —1, we have 

If we define A„ =^ J2iez+^ Ei^i^n{n G Z \ {0}), then we have the following Heisenberg 
relations 

[Am, A„] = m6m+n,oc. 

Proposition 1 The matrices E^ {(^ = 1, ( ^ 0, n E Z), A„ (n G Z, n ^ mod £) 
form a basis of 5le[s^^]. Therefore the representation tToo o L£ of sle is given by 

J2 EiX^""-' ^(A)^*(CA) : (C^ = 1, C 7^ 1), (2.7) 

ragZ 

A„ I— i> : ip-iip*^^ : (n G Z, 72 ^ mod i), Kg ^ id. 

Proposition 2 //we identify E^,Am with the elements o/gl^fs^^], t/ien we have 

s-Ei = Ei_,^„ .'"A„ = A„+m,. (2.8) 

Note that the centralizer of A^ in 0[(cxd) coincide with ie{sle) + CI. 
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Lemma 3 Let Qa{X) '= ipiX) ®^/'*(A) be the operator on JF^'" =^ Tj^®J^a- Then we have 



A=oo 



for Vn G Z. 



Proof. We use the expression 



For y4 = Ylii ^ij ■i'-ii'j- that satisfy ( |2.6| ), we have 

[A® 1 + 1® A$^V^-fc®^:+„^] 

k k 

k j k j 

= X] (^hk-n£ - aj+nl,k) 1p-j ® V'fc = 0, 

where we used ( PT^ ) in the last line. Q.E.D. 



(2.9) 



2.4 Equation for the SL^^'^-orbit of the vacuum 

Combining the results in the preceding subsection and Lemma |l], we have the represen- 
tation of sl^ on the space 

jrtor j,^ ^ defined 

^e°A = ('•^ ° TToo)*". Now we are in 

a position to define the following important operator: 

meZ 

The operator satisfy the following property. 
Lemma 4 We have 



J A=oo 



for any j G Z. 



Here the contour integral is understood symbolically, namely, just to extract the coefficient 
o/A-i : ^A"dA/(27rz) = 

Proof. We first note that the operator fi^''(A) is the product of fl^W = ip{\) ®'?/'*(A) and 
J2mez ^rn{X^) ® ^-m(A^) that commutc with each other. In view of lemma ^ and the fact 
that the latter operator is a series of A^, one can see the lemma above is valid. Q.E.D. 

Let SL^™ denote a group of invertible linear transformations on J^^"^ generated by the 
exponential action of the elements in sli ® R acting locally nilpotently. 
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def 



Corollary 1 Let r G J-'j^^ be in the SL^°^ -orbit of \a.Cj^ = vac^ vacj,. Then we have 



X^^n'XiX) (r ® r) rfA = /or Vj > 0. 



(2.10) 



A=oo 



tor 



Proof. By the preceding Lemma, it suffices to show the equations above for r = vac^ 



In fact, by the expressions (2.9) and (|2.4|), one can directly verify this. Q.E.D 



2.5 Bosonization 

Here we present a summary of the boson-fermion correspondence, that says the space JF4 
can be identified with the space of polynomials 

J^B = C[xi, X2, xa, . . . ] ® C[e"o, e"""], vacg = 1 G Jg. 
To state the correspondence precisely, we introduce the charged vacua 

vac^ n = • 

V'-n+i ■ ■ ■ V'-|vac_4 n > 

Lemma 5 /|^ There exists an unique linear isomorphism a : JF4 ^ jFg s'uc/i that 



a 



a(vacy) = va^^ crA^a"^ = <| 



m > 



mx-rn m < 

where we set vac^' =^ e"^° (n G Z). 
We introduce i/ie vertex operators by 

X(A) = exp(e(x,A))e"«A^^oexp(-e(4,A-i)), 

X*(A) = exp(-e(a;,A))e-^n-^-oexp(e(4,A-i)), 

X(A,/x) exp(e(a;,A)-e(a:,/i))A^-o^-9^oexp(-e(4,A-i)+e(4,/i-'; 

where dx stands for {dx^, \dx2, ^dx^, ■ ■ ■) and we set 

00 

e(:r,A) = ^x„A" (2.11) 

n=l 

Then we have the following correspondences of operators: 

a'^{X)a-^ = X(A), aip*{X)a^^ = X*(A), 
a(: V^(A)^*(/x) = ^(x(A,^) - l), 

where stands for the series Yl'^=o X^'^^^ji^. 
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2.6 Hirota bilinear equation arising from toroidal Lie algebras 

If we translate equation ( |2.10D into bosonic language, then it comes out a hierarchy of 
Hirota bilinear equations. So let us introduce the following bosonic counterpart of operator 

^TW = J2x{x)v^{x') ® x*{x)v^ux') = (fT® id)fi*7(A)((T-^ ® id). 

We define J-'^^ =^ JFg (g) jFg and vacg"' =^ vacg (g) vacg. To state next theorem, we define 
the elementary Schur polynomials Pn{x) {n = 0,1, . . .) by the generating series 

oo / oo 

^p„(a;)A" = exp ^x^A'' 

n=0 \fc=l 

here x stands for (xi, X2, X3, . . . ). 

Notation We shall use the notation y =^ (yi, y2£, ...),?/*== {y^, y^^, ■ ■ ■), y '= Vo, and 
y = {y, y) in the sequel of this paper. 

Theorem 1 Let t = rix.y) be in the SL\°^ -orbit 0/ vac^™ G J?^™. Then r satisfies the 
bilinear equation: 

00 

m,fc=0 

/or j > where a = (ai, 02, ...), 6 = 62^, . . . ) and k, are indeterminates and 

/^. = (a,Z).)^=^5^a„D.„, (6,D,)^=^^5^6„,D,„,. (2.12) 

n=l n=l 

Proof. By the bosonic realization, we shall identify (g) with the following space 
of polynomials 

C\e^^'\e^'"\x^,y^,y:!P {^ = 1,2, n = 1,2, . . .)]. 
To convert equation ( ^.10| ) into bosonic language, we shall introduce new variables. 





def 








def 




for 


n > 0, 




def 


-U^^ - 


y™^ J' 




def 


2 V^n^ ^ ' 


for 


n > 0, 


* 

yni 


def 


2 vyn^ 




yli 


def 


1/ *(1) , *(2)x 
Ipnt +yn\ ) 


for 


n > 1. 



Then we obtain 



X(A) ® X*(A) = exp I 2 ^ x„A" ) e'^oA^^o exp ( - ^ "TT^^^" 



1 / \ 1 '^dxn 

n=l I \ n=l 
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and 

J2 Kn(A') ® Vl„(A^ 



n=i / V „=i ^ 



Using a lemma by Billig ([0, Proposition 3, see also |T^), we obtain the Hirota bilin- 
ear equation. Here we made a reduction Dy*^ = 0, since, by the construction of the 
representation, r in the orbit does not depend on the variables j/*^. Q.E.D. 

Remark. 



1. By the construction of the representation, it is clear that our r = T{x,y) does not 
depend on the variables xe, X2£, .... 



2. If we put bi = = ■ ■ ■ = 0, then equation ( [^.121 ) is reduced to 



J2 Pm(2a)p„+,,+i(-D:,)e<"'^-V ■ r = (j > 0) 

m=0 

This is known as bilinear equation for the ^-reduced KP hierarchy 0. 

Example 2 Let £ = 2. From the coefficient of 03 for j = , we have a equation: 

[D^ - AD^Dt) r ■ r = 0, where x = Xi,t = X3. (2.13) 
This can be transformed into a non-linear equation 

2pxt = ^Pxxxx + Spl^ (2.14) 
where we set p = log r. Differentiating by x, we have the original KdV equation: 

ut = ^u^xx + ^uux, u = 2p^^. (2.15) 

Example 3 Let £ = 2. The coefficient of 62 in ( P.12| ) for j = gives an equation 



^Dy^Dl^ + ^Dy,D.,, - D,,Dy, ) T " T = 0. (2.16) 



Equation (|2.16|) times 6 is nothing but (|1.3|) with yo = y,xi = x, y2 = z,X3 = t. 



We shall use the following formulae of the logarithmic transformations: 



D/DyT-T D.,Dyr ■ r _ ^ 

2 ^^PxxPxy I ^Pxxxyi 2 ^Pxyi j 
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where p '= log r. These allows us to write equation ( p.l6| ) as a non-linear equation: 

_ 1 2 

2pa;j/2 '^Pxx Pxy ~l~ "^Pxxxy ~l~ g^?'*' 

Example 4 Let us derive Bogoyavlensky's 2 + 1-dimensional equation by our hierarchy 
of Hirota equations. We differentiate (p.l8|) by x to give: 



1 2 

'^Pxxy2 '^i,PxxPxy^ X ~l~ '^Pxxxxy ~l~ '^Pxyt (2- 19) 

'^PxxxPxy ~l~ '^PxxPxxy ~l~ "^Pxxxxy ~l~ g ^^^^Pxxxx ~l~ ^p^.^.^ (2.20) 

"^Pxxxxy ~^ '^PxxxPxy ~l~ ^PxxPxxy (2-21) 



We have used equation ( |2.15|) in the second line. Putting m = 2pa,a,, it reads 

Uy2 = ^U^xy + ^{d;^^Uy)Ux + UUy, (2.22) 

where we interpret d~^Uy as 2pxy. After the simple scale transformation dx = 2dx',dy2 = 
-dy'^,u = -4m', (12:2^ ) coincide with (1.22) in Chap. II in §. 



Example 5 Let i = 3. There is a non-trivial equation of degree 4 unique up to scaler, 
which is not in the 3- reduced KP hierarchy (the Boussinesq hierarchy), i.e. 

{ADx^Dy., - Dy,Dx, - Dy,DlDx2) r ■ r = 0. (2.23) 

For instance, this arises from the coefficient of 63 in ( |2.12| ) for j = 0. Introducing p = log r, 
we can rewrite ( p.23| ) as a non-linear equation: 



^Pxiys PyoX4 PyoxixiX2 ^Pyox\PxxX2 '2py^x2pxxxx — 0. (2.24) 



3 Lax formalism of Bogoyavlensky hierarchy 
3.1 Formal pseudo-differential operators 

The Lax formalism of the KP hierarchy is described in the language of formal pseudo- 
differential operators (PsDO for short) on a line. An affine coordinate x of this line is to 
be identified with the first variable Xi in the bosonization of free fermions in the previous 
section. Let dx denote the derivation d/dx. In an abstract setting, one can take an 
arbitrary ring with a derivation d (namely, a differential ring); we shall rather naively 
consider, e.g., the ring of formal power series of x. 

A formal PsDO is a formal linear combination, A = of integer powers of dx 

with coefficients a„ = an{x) that depend on x. The index n ranges over all integers with 
an upper bound N. The least upper bound is called the order of this PsDO. The first 
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non- vanishing coefficient is called the leading coefficient. If the leading coefficient is 
equal to 1, the PsDO is said to be monic. It is convenient to use the following notation: 

n>k 
n<k 

Addition and multiplication (or composition) of two PsDO's are defined as follows. 
Addition of two PsDO's is an obvious operation, namely, the term wise sum of the coef- 
ficients. Multiplication is defined by extrapolating the Leibniz rule 

to the case where n is negative. Here "o" stands for composition of two operators, and 
f^''^ the A:-th derivative d'^f/dx'^ of /. More explicitly, the product C — A o B oi two 
PsDO's A = Y,n «nC^x and B = Y,n ^n^^ IS givcu by 

m,n,k 

Note that the n-th order coefficient c„ = (C)„ is the sum of a finite number of terms. We 
shall frequently write AB rather than A o S if it does not cause confusion. 
Another important operation is the formal adjoint A i-^ A*: 

n n 

This is an anti-homomorphism, namely, for any pair ^4, B of PsDO's, 

{AB)* = B*A*. 

Any PsDO A = X]n<7v'^n^" with an invertible leading coefficient ojv has an inverse 
PsDO. In particular, a monic PsDO is invertible. 

3.2 KP hierarchy and its ^-reductions 

The standard Lax formalism of the KP hierarchy uses a monic first order PsDO (the Lax 
operator) of the form 

oo 
n=l 



(A) 



del 



>fe 



(A) 



def 



<fe 



iA)k 



def 
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The hierarchy is defined by the Lax equations 

dL 



dXr, 



with an infinite number of time variables x — (xi, 2:2, . . . ). The Zakharov-Shabat operators 
Bn are given by 

We identify xi with x, namely xi — x. This is consistent with the first Lax equation 

These Lax equations are associated with the hnear equations 

L*(A)=A*(A), ^i:^=B„*(A). 

Here ^'(A) (the "wave function") is understood to be a function of both x and A, ^&(A) = 
^(x, A), though we shall frequently omit writing x explicitly. The Lax equations ensure 
the existence of a non-vanishing solution of the above linear equations. Of particular 
importance is a solution (called the formal Baker- Akhiezer function) of the form 

def 



*(A)^^Ml + 5^m„A-" e^W, 



n=l 



where the first factor is generally a formal Laurent series with variable coefficients Wn 
Wn{x), and ^(A) is given by 



00 

def 



n=l 



The exponential factor is the same thing that we have encountered in the bosonization of 
free fermions. One can make sense of the action of PsDO's on ^(A) by simply extrapo- 
lating the derivation rule 

to negative powers of dx- 

Let us now introduce the Wilson-Sato operator 



def 

i + 

n=l 

This enables us to express the formal Baker- Akhiezer function as: 

*(A) - We^^^l (3.1) 
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The linear equations for \E'(A) can now be converted to the equations 



dW 

L = Wd,W~\ -- = B^W-Wd: (3.2) 

for W. The first equation shows the relation between the Lax operator L and the Wilson- 
Sato operator W (note that W is monic, hence invertible.) In particular, the Zakharov- 
Shabat operators can be written 



5„ = {Wd^^W'^) 



>o ■ 



If one substitutes the Zakharov-Shabat operators in the second equation for W by this 
expression, the outcome is the equation 

1^ = {wd:w-') ^^w- wo: = - {wd:w-') w. 

Note that this is a nonlinear system of evolution equations for the coefficients w„. This is a 
master equation of the KP hierarchy; the Lax equations can be derived from this nonlinear 
system via aforementioned relation ( p.2|) between L and W. It is this nonlinear system 
that Sato linearized on an infinite dimensional Grassmann manifold (Sato's universal 
Grassmannian) |2^, . 

The ^-reduced hierarchy (the £-KdV hierarchy) is defined by the constraint 

(LO<-i = 0. (3.3) 
The constraint means that be a differential operator, i.e.. 

This constraint is preserved under time evolutions, because the associated Lax equations 

dP 



dXr 



[(P"/')>o,P] 



for P become a closed system of unconstrained evolution equations for the coefficients 
M2, ... ,U£ of P. The time evolutions in xe, X2e, ■ ■ ■ are trivial, 

dP 

|^ = [P«,P] = 0, 

OXne 

and, in turn, give rise to an eigenvalue problem, 

P^(A) = A%(A). 

The other time evolutions can be interpreted as isospectral deformations of this eigenvalue 
problem. In the case where i = 2, this is exactly the well known characterization of the 
KdV hierarchy. 
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3.3 Heuristic Introduction of Bogoyavlensky hierarchy 

We here present a heuristic consideration that leads to a hierarchy of higher evolution 
equations for the Bogoyavlensky equation. 

The point of departure is the Lax representation 

dP dP 

^ = |pa„ + c.p| = p- + |c.p] 

mentioned in Introduction ( |1.4| ). Here P is the Lax operator + p of the KdV equation 
and C is a first order differential operator. This Lax equation is associated with the linear 
equations 

pvI/(A) = A2vl/(A), = {Pdy + C)^{X). 

Note that the second linear equation can be rewritten 

where 

dy 

Furthermore, if we assume the existence of the Wilson-Sato operator W with which ^^(A) 
takes the form 

^(A) = iye«W, ^(A) = a;A + --- , 

the time evolution of W is determined by the equation 

dW dW 
-— = P^ + CW 
oz oy 

= p[dy,w] + cw 

= {Pdy + C)W -Wdldy. 

We have used the identity PW = Wdl in the last two lines. 

If one compares the last equation with the equations in the KP hierarchy, one will 
soon notice that the role of 9^ is now played by didy. The (two-dimensional!) differential 
operator Pdy + C should be a counterpart of -B„'s. A natural generalization of this 
observation is to consider the higher order operators d'^dy, n = 1, 2, . . . . The associated 
evolution equations, with "time variables" ?/2(= z),y4,..., will accordingly take such a 
form as 



dW 



dy 



2n 



D2nW - Wdl'^dy 
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where is a two-dimensional differential operator in x and y. In fact, one can specify 
in more detail: This equation can be rewritten 

dW dW 

dy2n dy 

and this implies that — P^dy is a differential operator in x only, which we call C2n- 

D2n — P^dy + C2n- 

We are thus led to the evolution equations 

{p^dy + C2n)w-wdl''dy 



dy: 



2n 

dW 

= P^^ + C2nW. 

dy 

Just like the Zakharov-Shabat operators -B„ of the KP hierarchy, the differential operator 
C2n is uniquely determined by the evolution equation itself: 

f dW , dW ,\ f dW , 

\9y2n dy V 9y y 

Therefore the evolution equation of W can be rewritten 

dW f dW 1 » 



dy2n V 9y 



These equations, along with the aforementioned evolution equations in a;2n+i's, give a 
nonlinear system of evolution equations for the coefficients Wn of W . We call this system 
the Bogoyavlensky hierarchy or, more precisely, the Bogoyavlensky-KdV hierarchy . 

The spatial variable y may be identified with the zero-th time variable yo, because the 
corresponding operator Co is equal to zero so that the evolution equation for W reduces 
to 

dW _ dW 
dyo dy ' 

It is straightforward to write down the associated linear equations. They take two 
equivalent (but apparently different) forms. One expression is a direct consequence of the 
evolution equation of W: 



dy2n 

Another expression is the following: 



= (P"9, + C2n)*(A). 



dy2n 

where 



^(X'^dy+A2nMX), 



dy ■ 
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3.4 Lax equations of Bogoyavlensky hierarchy 

Let us derive Lax equations of the Bogoyavlensky hierarchy. To this end, we introduce 
the PsDO 

dy 

as the second Lax operator. This operator, just hke P = W o o W~^, is related to 
conjugation by the Wilson-Sato operator: 

dy-Q = W 0dy0W-\ 

The commutation relation [dy, 9^] = thereby implies the commutation relation 

[Q-dy,P]=0 

or, equivalently, 

dP 

Note that this relation itself takes the form of a Lax equation. 

These operators have already appeared in the aforementioned differential operators 
and A2n- Can be indeed written 

Moreover, using the Lax-type equation 

Qpn 



dy 

(which is a consequence of the commutation relations of P and Q) , one can confirm that 
can be written 

A2n = -(<5P")>0. 

Having these operators, we can now write down Lax equations of the Bogoyavlensky 
hierarchy: 

Proposition 3 The operators P and Q satisfy the equations 

[<9x2„+i — -B2n+1, -P] = 0, [dx2n+l ~ B2n+l, Q — dy] = 0, 
[dy.n - P'^dy - C2n, P] = 0, [dy,„ - P^dy - C2n, Q - Oy] ^ 0, 



or, equivalently, 



dP 



dx 



2n+l 



dQ dB2„+i 



dx2n+i dy 
dP „^ dP 



[-B2n+1, P], 

+ [-B2n+1, Q]-, 



P"l^ + [C2n,P], 



dy2n dy 



dy2n dy dy 
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Proof. Let us prove the equations for the ?/2n-derivatives; the proof of the other equations 
are parallel and even simpler. First, differentiating P = Wd^W'^ gives 



dP 



dlw-^-wdlw~^^w-^ 



2n 



dy2n 
dW 



dy: 



2n 



dy 



w-\p 



2ri 



We can now use the evolution equations 

dW 



pnqw + C2nW 



2n 



and the commutation relation [Q, P] = dP/dy, and find that 

dP 



dy 



2n 



[P^Q + C2n,P] 
P-^ + [C2n,P]. 

dy 



Similarly, from the identity Q — dy 



-WdyW-\ 



dQ 



dy: 



211 



dW 



W-\Q-dy 

dy2n 

= [P'^Q + C^n.Q-dy] 

The first two terms in the last line cancel each other, again because of the commutation 
relation of Q and P. This completes the proof. Q.E.D. 



3.5 Formulation of £-Bogoyavlensky hierarchy 

We are now in a position to formulate the £-Bogoyavlensky hierarchy for a general value 
of £ > 2. This hierarchy is an extension of the ^-reduced KP hierarchy with additional 
independent variables y = {y,y), where y = yo and y = {ye, y2i, ■ ■ ■)■ The coefficients Wn 
of the Wilson-Sato operator 

oo 
n=l 

are understood to be functions of these variables also, Wn = Wn{x,y). Assume that W 
satisfy the £- reduced constraints, cf.( p.3|) : 

{wdiw-')^_^ = o, 



20 



and define the differential operator 

P ''^ wdiw-\ 

• Tfie most fundamental part of this hierarchy is the evolution equations 

dW 



dXr 



- -{wd:w-')^_^w, 



dyne V dy y 

of W. These equations can also be written 
dW 

dyne oy 
where and Cne are differential operators defined by 



dy 



>o 



• The evolution equations of W in ?/„^'s can also be rewritten 

dyne dy 

where 

d f dP^ 
A-ni = — (Q-P")>0 = Cne 7{ — • 

dy 

• Lax equations of this hierarchy are described by P and the following PsDO 

dy 

They satisfy the commutation relation 

dP 

Bn and Cne can now be expressed as 

Bn = (P"/0>o, Cne = -(P"g)>o- 
21 



P and Q satisfy the Lax equations: 

dP 



dXn 

dQ OB 



dxn dy 
dP „„dP 



- + [Bn,Q], 



dyne dy 
dQ _ dQ dCni , . 
dyni dy dy 

The formal Baker- Akhiezer function 

vl/(A) •^=?:Ve«W (3i 

satisfies the hnear equations 

d^{X) 



5y 



^^=Bn^{\), ^^ = [p-dy + Cni)^{X). (3.9) 



The last equation can also be written 



^l^=[\-% + An,)^{\). (3.10) 



Example 6 Let us consider the case £ = 2. The Lax operator has the form P = d^+u 



where u = —2wix- We can write down C2 = A2 + ^ as follows: 



C2 = Cidx + C2, Ci = -Wiy, C2 = WiyWi - W2y ~ 2Wir,y. (3-11) 

The evolution equation with respect to y2 reads 

dP dP 

^ = P^ + [C2,P]. (3.12) 
dy2 dy 

One can calculate the right-hand side of (|3.12|) as follows: 

P^ + [C2,P] (3.13) 

= {2Uxy - Cij:x - 2C2x)dx + Uxxy + UUy + CiUj: - C2xx, (3-14) 

where the coefficients of d"^ cancel trivially because u = —2wix- 

Here we use the following relation which is a direct consequence of the 2-reduced 
condition 

wixx + 2w2x - 2wiWix = 0. (3.15) 
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Using this, we now have C2x = —^Wi^xy This finally leads to 

^— = -^Uxxy + vux + uuy where v =^ —Wiy. (3.16) 

Caution! The expression W2x, for example, does not mean d^w but dx'W2- 

Example 7 We shall demonstrate the third order flow in case of £ = 3. The Lax 

operator P = + Udx + V is given as: 

U = -SWix, V = -SWixx - 3W2x + ^WiWix- (3.17) 

Operator C3 = Cid^ + C2(9^ + C3 is given explicitly as: 

Ci = -Wiy, C2 = WiyWi - W2y - SWixy, (3.18) 
C3 = -W3y - WiyWj + Wi{w2y + 3Wixy) +Wiy{w2 + 5Wix) ~3Wixxy - 3W2xy-{^-'^Q) 

Now we recall the 3-reduced condition. In particular, the coefficient of d^^ in Wd^W'^ 
should vanish identically. This condition is written as: 

U)3x = -^Wlxxx - W2xx + + ^i^Wa + WiWixx + ^1^2^ " ^Ix^^ (3.20) 

Using ( |3.171 ),( pTT8D , and ( p.l9| ), it is simple to see that the right-hand side of evolution 
equation, 

dy^P = P^+\C,,Pl (3.21) 

is of order at most two. To verify that the right-hand side is actually of flrst order, we 
eliminate W3 by using equation ( |3.20|) . Thus, it should be emphasized that the 3-reduced 
condition is crucial to introduce flow y-^ by ( |3.21D . 

After eliminating w^, we have the following explicit forms of evolution equations: 

dy^U = 6WiyW2xx - Wixxxxy + QWlxyW2x " '^W2xxxy " 12u^lyU^L - 15WixyWiWix 
+6WixW2xy - 9WiyWiWixx + l2WixyWixx + QWixxyWix + 2WiWixxxy 
+5WiyWixxx + 3WixxW2y, 

(^ysY = -^WixyWl^ - ISWixyWiWixx + ^WixxxW2y + ^WixW2xxy - ^wl^W2y 
-15WiyWixWixx + QWixyW2xx + QWixxyW2x " l2WixxyWlWix 
+3WiyW2xxx - QWiyWiWixxx + 2WiyWixxxx + 7WixyWixxx + ^WixxyWixx 
+3WixxxyWlx + Q'WixxW2xy " 3W2yWiWixx + l2WiyWiwl^ + SWiyWjwixx 
-9W2xWiWixy + 9W2xW2xy + 3W2yW2xx " 9'W2xWiyWix - 3WiyWiW2xx 
QWiWixW2xy ~l" QuJiWixWixy U]2xxxxy 3'U^lxxxxx2/ ~l" W\W\xxxxy 
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3.6 Explicit forms of operators 

We give a supplementary discussion on the forms of differential operators An- The dif- 
ferential operators Ani and Cne are written in terms of W (|3.7| ) , (|3^ ) and hence in the 
coefficients Wj's of W. In particular, we have 

Here we shall recursively determine the coefficients a^p\ If we cast ^E'(A) = We^^^^ into 
linear equation ( p.lO| ), then we obtain 

dyn dy 
Taking the polynomial part of the equation above with respect to A, we obtain 

n / n \ n 

i=l \j=0 / j=l 

Equating coefficients of A'^"^ in the equation above, we obtain a^^^ = —dyWi. Moreover, if 
^ , . . . , are already known, then a,-"'' can be uniquely determined by the coefficients 

of A"~*. In particular, we have Og"'' = widyWi — dyW2- Note that af^ is a differential 
polynomial of wi, . . . ,Wi with respect to dx and dy. 

Example 8 For any n we have of'"' = —wiy, = wiyWi — W2y We shall give the 
formulae for n < 5: 

ai'' = -WSy - WlyWl + WlW2y + WlyW2 + 2WlyWlr,, 
4'^'' = WlW2y - WlyWl + SWlyWlr, - WSy + WlyW2, 

(4) 2 S 

04 = -Wiy - W^W2y + WlyWl - hWlxWlyWl + WlW^y " 2W2WlyWl 

+W2W2y + 2WixW2y + WiyW^ + 'iWiyW2x + ^WiyWi^x, 

aP = -WlyWf + WlW2y - WSy + AWlyWlx + WlyW2, 

af^ = QwiyWixx + wiywf - wfw2y + WiWsy - IwixWiyWi - 2W2WlyWl (3.22) 

+3WixW2y + W2W2y - Wiy + AWiyW2x + WiyWs, 

(5) 2 
Og = -W^y + WiyWi + WsW2y + 3WixxW2y + 3W2xW2y + ^WiyWixxx + 'iW2WiyW^ 

-2w2WiW2y - 6W2WiyWix + 9WixWiyWl - 5WixWiW2y - 2w3WiyWi - 9WixxWiyWi 

-7W2xWlyWl + 6WlyW2xX + 4:WlyW3x + ^2^3?/ - WlyWj + 2WlxW3y " SWlyWl^ 
-WlyWl + Wlw2y - Wjwsy + WlWiy. (3.23) 
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4 Reproducing bilinear equations in Lax formalism 

4.1 Dual formal Baker- Akhiezer function 

We define the dual ^'*(A) of tlie formal Baker- Akhiezer function ^(A) as follows: 

Proposition 4 \E'*(A) satisfies the linear equations 

P*^(A) = A^^^*(A), g*^(A) = -^^^A^, 
^^^^ = -b:^*{x), = (P*"a, - A^,) **(A). 



The last equation can also be written 



dyne 

Proof. The first and second equations are immediate from the operator relations 

oy oy 
The third equation can be derived as follows: 

dW*-'^ 1 dW* 



Similarly, 



dxn dx. 



dW*-^ 1 dW* 



. w*-'-—w 
oynt oyne 

dW 

a ^ 

oynt 

dy 

which implies the fourth equation. The last equation can be readily derived from the 
fourth equation. Q.E.D. 
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4.2 Residue formula of product of PsDO's 

The following lemma 0] is a clue for deriving the bilinear equation that ^E'(A) and 
^*(A) satisfy. 

Lemma 6 For any pair A, B of PsDO's and non-negative integer j, 

2^ i oo ^^'"''^ ^ {-iy^\BA*)-,-,. 

In particular, if A and B are PsDO's of zero-th order of the form A = + 0{d~^) and 
B = h, + 0{d-^), 

BA* = boao ^ (Vj > 0) / {diAe''^) (Ee"^^) dX = 0. 

J A=oo 

Proof. Let A and B be written A = a„9" and B = J2n ^nd^- Since 



m k>0 

the contour integral can be calculated as: 



k—m—n—l=j 

On the other hand, 



Tri /a=oo 



2 

27ri 



m,n 
m,n,k 

m,n,k ^ 
J fc— m— n— l=j ^ ^ 



Therefore the coefficient of d^^ ^ coincides with the countour integral. The lemma is thus 
proven. Q.E.D. 
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4.3 Bilinearization of linear equations 

Let us slightly change the notation, namely, we make explicit the functional dependence 
on X and y, e.g., '^{x,y, A), \E'*(x,|/, A) and ^{x, A) stand for \E'(A), \E'*(A) and ^(A) in the 
previous notation. 

We now apply Lemma ^ to the case where A = W and B = W*'-^. This eventually 
leads to a bilinear equation for \E'(a;,y, A) and '^*{x,y, A) as follows: 

1. By the lemma applied to A = W and B = W*~^ gives the bilinear identities 

^_ {di^{x,y,X)^^*{x,y,X)dX = 
for j > (note that X2A^ + x-^X^ + ■ ■ ■ in e^'^'^^'^^ cancel out). Thus, we have 

jf (^m{x,y,X)')-^*{x,y,X)dX = 

for any differential operator R = J2n>o''^ri{x,y)d!^ in x. 

2. Iteration of the evolution equation of \E'(a:,y, A) gives rise to higher order equations 
of the form 

d:id:^^ ■■■^ix,y,X)= 5a„a„...vl/(x,|/, A) 
for «!, • • • > 0, -Bai,a2,... being a differential operator in x. 

3. Combining these equations with the last bilinear identity, we obtain the bilinear 
equations 

£_ (^."^C ■ ■ ■ ^(^'^' A))m/*(x,|/, A)rfA = 
for «!, a2, . . . > 0. 

4. These bilinear equations can be cast into a generating function: Introduce an infinite 
number of new variables a = (oi, 02, . . . ) , and sum up the bilinear equations over 
«!, ^2, . . . > with the weight a"^a2^ ■ ■ ■ / ai!a2! ■ ■ ■ • The outcome is a single bilinear 
equation of the form 

j) <i/{x + a,y,X)'^*{x,y,X)dX = 0. 

J A=oo 

5. Finally, rewriting a; + a to ac', we obtain the bilinear identity 

S ^(x', y, A)^*(x, y, X)dX = 0. 

J A=oo 
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Remark. Note that the status of xg, X2e, . . . are somewhat different from others: Since 
Wn do not depend on these variables, they can appear in e^^*^^'^^ only. The net effect of 
these variables is thereby insertion of the exponential factor 

oo 

exp(^^«^-a;„^)A"^) 

n=l 

in the contour integral of the bilinear equations. Therefore the last bilinear equation can 
be further reduced to the bilinear equations 

J A=oo 

for j > 0. 

This is not the end of the story. What we have done is simply to bilinearize the KP-hke 
part of the £-Bogoyavlensky hierarchy. The other evolution equations i/nes still remain to 
be bilinearized. 

Bilinearization of the evolution equations in ynes can be achieved by the following 
steps, which are almost parallel to the previous calculations: 

1. We rewrite the evolution equations as 

- X''%)^{x,y, A) = Ar,t^{x,y, A) 

and again do iteration. This yields the higher order equations 

{dy,-X%f^{dy,,-X'%f^---^{x,y,X) 
= Ap^^p^^,„^{x,y,X) 

for Pi, P2, • • • > 0, A^i,/32,.-- being a differential operator in x. 

2. Combining these equations with the bilinear equations that we have derived above, 
we now obtain the bilinear equations 

{{dy, - >^'dyf'{dy,,, - X'%r^ ■ ■ ■ ^{x',y, A)) 
X <if*{x,y,X)dX = Q. 

3. Introduce new variables b = {bi, 62^, . . . ) and sum up the bilinear identities over 
1^1,^2, ■■ ■ > with the weight h^f^h^l ■ ■ ■ / Pi\(32 \ ■ ■ ■ ■ This sum contains power series 
of dy^^ — X^^dy, which is essentially an exponential: 

/3>0 n=l n=l 
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The action of this exponential operator can be easily understood by the obvious 
identity 

oo oo 

exp(5^6„,(9,„, - X-%))f{y,y) = f{y - ^ 6„,A-^ j/ + ft) 

n=l n=l 

that holds true for any function f{y,y). We thus obtain the bihnear identity 

^ (x', y - 7j{b, X),y + b,x) X)dX = 0, 

where for indeterminates b — {bi,b2e, ■ ■ ■) we set 

oo 

r,{b,X)'^Y.bnA'''. (4.1) 

n=l 

4. These calculations can be repeated, now starting with the linear equations for 
'^*{x,y, X). This yields a bilinear equation with the variables y = {y,y) in ^* 
being shifted as 

{y,y) — ' {y-v{c,x),y + cj 

where c = (q, C2e, ■ ■ ■ ) are newly introduced variables. 
We thus eventually find the following result: 

Theorem 2 The formal Baker- Akhiezer functions ^{x,y,X) and ^*{x,y,X) satisfy the 
bilinear equation 

jf ^(x',y-ri{b,X),y + b,X^ 

X '^*(^x,y-ri{c,X),y + c,X^dX = 0. (4.2) 

Here x, x' , y, b and c are understood to be independent variables. 
Remark. 

1. Conversely, one can reproduce the evolution equations of W etc. from these bilinear 
equations. 

2. The preceding remark on the status of X£, X2e, ■ ■ ■ also applies to this bilinear equa- 
tion. Thus the bilinear equation holds true if any nonnegative power of is inserted. 
Accordingly, by taking a linear combination of those equations, one eventually ob- 
tains the slightly generalized (but actually equivalent) form 

jf_ f{X')^(x',y-ri{b,X),y + b,x) 

X '^*(x,y -r]{c,X),y + c,X^dX^O. 
of the bihnear equation that holds for any power series /(A^) = ^„>o/nA"^. 
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4.4 r function and Hirota bilinear equations 

We can define the r function in tlie same way as tlie case of the KP hierarchy. Namely, 
the r function is a function r = r(x.y) that satisfies the equations 

A) ^ ^^^;^e«-., (4.3) 
T{x,y) 



where 

1 1 1 
A' 2A2''" 



e(A) = 



Note that this definition allows the indeterminacy 

T{x,y) — > f{y)T{x,y) 

with f{y) being an arbitrary function of y and y only. Equation ( [4.2|) now turns into the 
equation 

T[x' - e{X),y - Tjih, X),y + 6) r (ar + e{X),y - r]{c, X),y + c) 

X ^— ^ ^ ^±dX = 0. 

r(x',y -r]{b,X),y + hj T\x,y - r]{c, X),y + cj 

for the T function. Notice here that the two factors in the denominator are power series of 
A^. According to the second remark of the last theorem, one can insert any power series 
/(A^) in the foregoing equation. If we choose 

/(A^) = T(x',y- r](b, X),y + b^T(x,y - r]{c, X),y + cj, 

the two factors in the denominator are cancel out. The outcome is the bilinear equation 

exp^^(:r' — X, A)j 

X t(x' - e{X),y- r]{b, X),if + b^T(x + e(A), y - r]{c, X),y + cjdX = 0. 

for the T function. 

To rewrite this bilinear equation into the Hirota form, we replace 
and choose c to be equal to —6. The bilinear equation can be thereby converted into a 
Hirota form: 



/ 

J X=oo 



I exp (2 ^ a^X^^ exp (- ^ ^ A"") exp (- ^ h^,D. 

•^^=°° n=l n=l n=l 

X e<"'^->+<'''^« y) ■ t{x, y)dX = 0. 
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The exponential including A can be expanded into a series of Schur functions: 



exp(2 5^a„A") = 5^p„,(2a)A", 

n=l n=0 
oo p. oo 



n=l n=0 

oo 



exp(-5^6„,A"^Z},) = 5^p„(-6D,)A"^ 



71=1 n=0 

By plugging these identities into the last bilinear equations and do contour integrals, we 
obtain bilinear equations of the form 

oo 
m,k=0 

X e<'''^->+<^'^*>r(a;,|/) ■ r(a;,j/) = 0. 

These bilinear equations almost coincide with those that have been derived from the 
representation theory of toroidal algebras. Although they do not agree, the discrepancy 
is rather superficial: If one shifts a to a + K,b/2 and move the extra exponential factor 

exp 6„£A"^ j to the third exponential including Dy, the outcome exactly reproduces 

the previous result (cf. Remark in the preceding subsection). 



5 Special solutions of hierarchy 

5.1 Construction of the solutions of the Wronskian type 

We shall apply the method developed by Date in (||^, ^, [T^) to construct special solutions 
of £-Bogoyavlensky hierarchy. The method corresponds to consider Baker-Akhiezer func- 
tions on rational singular curves with nodes. Such algebro-geometric backgrounds are 



given by Krichever in (see also of Appendix 1 for an exposition) and Manin in 
2^ Chap. III. One can also consult instructive lectures given by Mumford in pSf (Illb, 



section 5) that treat the subject. 

As the data for the solution constructed below, let us consider the following: 

• ai = ai{y) , Ci = Ci{y) {i = 1, . . . , N) are arbitrary (local) solution of the equations 

af— = 0, aTj- = ^ n = l,2..., 5.1 

oyni oy dyne oy 

such that we have ^ aj for i ^ j, 

• (i {i = 1, . . . , N) are i-th roots of unity not equal to 1. 

• Tj (i = 1, . . . , A^) are arbitrary complex numbers. 
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We put l3i{y) =^ (iai{y). Then we have the equations 

dyne oy 



Note that equation (5J) have constant solutions, namely, for distinct complex numbers 
ai {i = 1, . . . , N) and arbitrary complex numbers Cj {i = 1, . . . ,N),ai{y) = ai, Ci{y) = Ci 
trivially satisfy equations ( ^.1|) . 

Let us construct a special solution of £-Bogoyavlensky hierarchy, which we shall seek 
in the form 

^(A, r) = (A^ + w7iA^-^ + --- + wn) e^^(^), (5.2) 

with Wn = Wn{x,y) being unknown functions. Here we set '= ^{^^ X) + ry + rri{y, A) 
for any constant r G C (cf. (|2ll|) ,(^). 
We define the N x N matrix 

S;v=(/(°\ /«,..., /(^-^)), f^'^dlf 
where / denote the iV-column vector . . . , fjy) and 

f^^^ ^irM) + c,e^rM). (5.3) 

Note that, for any constant solution {ai, q} indicated above, we clearly have det En ^ 0. 

Theorem 3 Let {aj,Cj,Cj,rj {i = 1,... ,A^)} be a given set of data described above. 
Suppose we have det Hat ^ 0. Then the condition 

^(A,r,)+c,^(a„r,) =0 {t = 1, . . . , N) (5.4) 



uniquely determine the function \l/(A,r) of the form in ( \5. Sj ), and '^{X) = \l'(A,r) solves 



equations ( \3.1C\) . Explicitly, we have the following expressions 



det(/W,...,/(^-^),/W,/(^'+^),...,/(^-^)) 

det(/(o),... ,/(^-i)) ^ i,---,^vj. 

In particular, the corresponding r function is given by the formula 

T = det Sat. (5.6) 

Proof. Using the obvious differential equations 

the condition (|5.4|) can be written in the following form 

WnI = (5.7) 
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where by Wn we denote the differential operator 

Wj^ = + w,d^~' + --- + wj^. (5.8) 
Since detS^v ^ 0, we can solve equation ( pT7| ) by Cramer's formula, to obtain expression 

To show the function \E'(A,r) satisfy (|3.9D, we shall show that the PsDO defined by 

W = WnO-^ solves equations Q and 

It suffices to show equation (|3.4| ) for some differential operator Bn in x, since differen- 
tial operator _B„ is uniquely determined by the evolution equation itself. Differentiating 
with respect to x„ and using dx„f = i9"/, we have 

OWn \ 

+ WnO: / = 0. (5.9) 



There exist differential operators and R in x such that 



+ WNd^ = BnWN + R, ord(i?) < AT - 1. (5.10) 

OXyi Ox 

/^From ( ^.7|) and ( ^.9|) we have 

Rf = ii = l,...,N). 

This implies R = since detS^r ^ 0. Hence by multiplying d~'^ from the right to 
equation ( |5.1CI| ) with i? = 0, we obtain ( |3.4| ). 

Now we recall a^ = Since the dependence of ft on Xni comes only through the 
overall factor exp [a^^Xne) = exp , the expression ( |5.5| ) imphes 

dW 

^ = (n = l,2...). 

Hence from ( |3.4|) Wd^^W~^ is a differential operator, namely we have 

Wd'^^W-^ = WnO^W^^ = P", (5.11) 

where P denotes the differential operator Wd^W'^. 

Next we prove ( p.5| ), the evolution equation of W with respect to Using ( p.l| ) and 
( ^.3|) one has the equations 



dy.J={d:%-^)f for n=l,2,..., 
where $ diag(a,«), . . . , a,«))- 

In view of these equations, the differentiation of (|5.7|) with respect to yields 



+ WNdfdy - Wn<I>^ / = 0. (5.12) 
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We use (pi 



) to rewrite the operator W^dyd^^ as follows 



dy 



N 



WNdydf = WNdfW^^WNdy = P'^WNdy = P" 
By virtue of (|5.13| ), Wn^ = ^Wn and (|57|), we can rewrite ( ^.12|) as 

/ = o. 



+ dyOWr, 



(5.13) 



dWN _pndWr, 



dyn 



dy 



(5.14) 



Note that the operator in the left hand side of ( ^.14| ) is an ordinary differential opera- 
tor, which does not include dy. So we can apply exactly the same argument above to get 
the unique differential operator Cne in x such that 

oyni oy 

This completes the proof of ( p.5| ). To see the formula ( |5.6|) is valid, we only have to notice 



(5.15) 



the following identity arising from ( |4.3|) and ( p.ll) : 

wi = -(9a; logr. 

Q.E.D. 



5.2 A^-soliton solutions 

li ai, Ci{i = 1, N) in the preceding subsection are constants, namely do not depend 
on y, then the solution of Wronskian type has clear representation-theoretical meaning. 
In fact, we shall see that, up to an irrelevant factor, =^ detS^r coincide with the 
following N-soliton solution 

4 = (1 + ar,X{aN,MVmA(^%)) ■■■{! + aiX{ai, (3i)V^M)) va4°^ (5.16) 

where are constants we shall specify later. To be rigorous, we have to assume |aAr| > 
■■ ■ > I ail so that the series converges. The fact that is a solution can be readily 
observed in view of the next lemma (cf. Lemma 

Lemma 7 Let = 1, C 7^ 1 o,'>T'd m E Z. Then on the representation {J-'q'^, vr^°g) we have 
^£;^rA-" ^ ^(X(A,CA) - l)Vm{X'). (5.17) 

Proof. Let us put 

r(A) = Eit- ■ X- + j^K:^{X'). (5.18) 
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Then we have 

[(j)l,r{X)] = mA^^r(A) VfcGZ (5.19) 

and 

[Afc,r(A)] = (1 - c'')A'=r(A) Wkez. (5.20) 

( p.l9| ) is a easy consequence of ( |2.8| ), and ( p.2(]| ) follows from the following relation, which 
can be directly verified 

[Afc, ] = (1 - C'^) {Ei^kt"" + log VA; G Z (5.21) 

where 6m,n = 1 if m = n mod i and otherwise. 

By virtue of a standard lemma on vertex operators (see e.g. |TB| Lemma 14.5), we 



have the conclusion that r(A) is realized on JF^™ by 

c^,mX(A, CA)Kn(A^), Cf^m is a constant. (5.22) 

To determine the constant c^,m, we notice that EQf^vac^^^ = 0. To see this, we only need 
E'gvacg = 0, which is obvious from (|2.7|), (p.5[), and the description of the representation 



^0 

in lemma |l[ Now it is immediate to see c^^m = (1 ~ C)~^) ^^^id the proof completes. (J. 
If we define the following matrices 

E, = diag(e«("i),... ,e«("-)), C = diag(ci, . . . , c^), = {a{-\j=,„„,^. 

then En = E^V^ + C EoMoi and hence we have 

detS;v = detE^detV>det(l + Dr), D = CEa,E^\ T = Va,Vp\ (5.23) 

For any matrix X = {xij)ij=i^,,,^N and any subset J of {1,2, . . . , A^}, we define the 
principal submatrix Xj = {xjj/)jj'(zj. We shall use the well-known formula 

det(l + X)= det Xj, J runs over all the subsets of {1,2,..., N}. (5.24) 

JC{l,2,...,Af} 

Lemma 8 Let T =^ VaV^^ = {'yij)i.jei- Then we have 

Combining the formulas ( [5. 23] ), ( p. 24] ), and the above lemma (cf. |^3|), we have 



Jci VieJ / ij'e. 



J 



where A(^) denotes Vandermonde's determinant det = ni>j(A ~ Pj)- Here we also 
note that does not depend on ?/*^ from the construction. So if we set Oj = Cj'jjj, we 
can see that rjf is equal to times an irrelevant factor e^'=i^™''^'^'^A(;9). 



35 



6 Conclusion 



Inspired by the work of Billig and lohara et al., we have introduced the £-Bogoyavlensky 
hierarchy and clarified its group-theoretic properties, Lax formahsm, and special solutions. 
The group-theoretic characterization of the new hierarchy is based on the 2-toroidal Lie 
algebra slf^ . We have constructed a representation of this Lie algebra, and derived a 
system of Hirota bilinear equations for the r function. The point of departure of the 
Lax formalism is the observation that the lowest two members of the Hirota bilinear 
equations of Billig and lohara et al. are equivalent to Bogoyavlensky's 2-1-1 dimensional 
extension of the KdV equation. Bearing this in mind, we have constructed a system of 
Lax equations, from which we have been able to reproduce the Hirota bilinear equations. 
The Lax formalism has also turned out to be useful for understanding special solutions. 

Several interesting problems remain open. Firstly, the present construction should 
be extended to the toroidal Lie algebras associated with simple Lie algebras other than 



sli- This will give an extension of the Drinfeld-Sokolov hierarchies in the spirit 
of Bogoyavlensky. The same problem can also be raised to another important family 
of soliton equations, namely those represented by the nonlinear Schrodinger equation. 
Strachan's work p9| is remarkable in this respect. He presented a 2 -|- 1-dimensional 
extension and an associated hierarchy of higher equations for soliton equations of this 
type. Lastly, we would like to mention the problem of constructing special solutions of 
the algebro-geometric ("finite-band") type. A precursor of this problem is Cherednik's 
work 0, in which he presented a construction of algebro-geometric solutions to the self- 
dual Yang-Mills equation. Since our hierarchy is closely related to the self-dual Yang-Mills 
equation, a similar construction of special solutions will be possible. 



7 Appendix 

We give a list of the Hirota equation of low degree contained in the hierarchy for i = 2. 
Since we have P{D)f ■ / = for any function / and polynomial P such that P{—D) = 
—P{D), we shall list below only the even polynomials in D. We also drop the terms 
including D^^, D^^, .... We shall use the abbreviated notation Dn and Qn for Dx„ and 
Dy^ respectively. 
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degree 4 




degree 5 


('n^ 1 on n^n j.oAn'in i9nnn 


( Ji' -J- 90 n- n ^ 9 1 n n foo n-^ ^ inn d 


(Z^i^ - lOD^Z^a + Qo + (lODt - AOD3) Q2 


degree 6 




nt) '?9 _i_ 4 nii n„ 


nt* snn^ 9n n-^n j^aaati n 


degree 7 


^'n'''-L79^n j-vnn^n j-osn n^n j-^^nAn n^^ /n 
- (420^2^3 + SlZ^i^ + 504D5) Q2 - (420^3 ^ 840/^3) 


('in^ 1 Rs n n^ 1 /icnn /n nionn osnnf'i/n 


("i n 1 fis n n^ 1 /i cn n ^ /n j_ ^'osn n^ n os n^ fi7o n ^ /n 


{dDi - IUOD7 - rODfDa + SeODpi) Qo - (126D^ - 2OI6D5) Q2 


( -1. 790 n™ -I- 7n n* n„ -1- 9sn n^ n, -1- "^rtA n^ n^^ r)„ 

- (840^2^3 + 42Di5 + 10O8D5) Q2 


{Dl + 72OD7 + YODfDs + 280D^Di + 504L)5Di^) Qq - 50A0DiQq 


\d\ + 72OD7 + 1QD\D^ + 280^^1 + hQAD^D'i) Ql - 5040DiQ^ 
- (2520^2^3 + 12QDI + 3024D5) QIQ2 + (2520^3 + 504OD3) QoQl 
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